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Abstract 

We discuss the limitations of the standard coUinear approach. The 
kinematical approximations necessary to derive the colUnear factor- 
ization are insufficient for the description of the exclusive final states. 
We argue that for a proper treatment of the final states one needs to 
use fully unintegrated parton correlation functions. We introduce the 
gauge invariant definitions of these objects and the factorization theo- 
rem for one jet production in deep inelastic scattering. 



1 Introduction: factorization and integrated parton distribution functions 

The leading twist formalism based on the collinear factorization [1] has been proved widely suc- 
cessful in a variety of processes in QCD. The collinear factorization allows for the separation 
of the short and long distance contributions in the processes which involve hadrons. The short 
distance part is the partonic cross section, which can be systematically calculated order by order 
in the strong coupling within the perturbation theory. The long distance parts are parton distri- 
bution and fragmentation functions which contain all the non-perturbative information. Thanks 
to the factorization theorem these quantities are universal and once determined in one process 
can be taken over to the other process to determine the corresponding cross section. These stan- 
dard parton distribution functions are often referred as the integrated parton distribution functions 
(PDFs) since they depend only on the fraction of the longitudinal momentum x and the hard scale 
jjL. Problems arise when one tries to address more exclusive processes with final states. Then, 
the details of the kinematics are essential and standard kinematic approximation can be insuffi- 
cient for these purposes. This is where the formalism with the parton distributions unintegrated 
over the other components of the momenta is necessary. One also needs to prove more general 
factorization theorems which incorporate these unintegrated parton distribution/fragmentation 
functions. 

To start with, let us first recall the standard definition of the integrated quark parton density: 



Here, %l) is the quark field and \P) is the proton state in which the operator above is evaluated. 
Subscript R means that we are considering renormalized operator and ji is the renormalization 
scale. The object inserted between the quark fields is the Wilson fine 
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in the lightlike direction n = (0, 1, 0^)- This ensures the gauge invariance and can be shown to 
arise from graphs with arbitrarily many gluon exchanges which appropriately factorize into this 
path-ordered exponential. 

To illustrate some issues related with the kinematical approximations used in the standard 
coUinear approach let us consider the parton model of the deep inelastic scattering. The lowest 
order graph is pictured in Fig. [l] and we can write down formally its contribution to the hadronic 
tensor as ^ 



47r J (27r)^ 



with the sum over j as a sum over the quark flavors. 




Fig. 1 : Parton model in deep inelastic scattering. 

The kinematic approximations needed to obtain the factorized form set the quarks on-shell 
and massless 
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These approximations give the familiar expression for the hadronic tensor 



4tt 



dk (fh 
^2^ 



(4) 



(5) 



The expression in the curly brackets corresponds to the parton density and the expression with the 
trace is the partonic cross section. Note, that the approximations are good when the momentum 
components kr and k~ are small with respect to k~^. This is sufficient when the momenta are 
integrated over in the inclusive cross section but might lead to large errors when we look at 
some more exclusive processes. Indeed the parton model approximation makes the following 
replacement 

6^^\k + q-l)^ 6{k+ + q+)5{q~ - 1-)6^'^\It) ■ 



These two delta functions can lead to similar results when integrated over all momenta (inclusive 
processes) but can lead to very different results in the exclusive cases. If we allow for the invariant 



mass of the outgoing parton P = we see that the parton model value = xP^ becomes 
significantly modified 

k+ = a;P+ + ^ 
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As shown in [2] on the example of the cc production in DIS there are huge differences be- 
tween standard parton model, unintegrated partons and exact kinematics. In particular the pT 
distributions of the produced charmed quarks depend heavily whether integrated or unintegrated 
formalism is used. The general conclusion [2] is that conventional formalism with integrated 
parton densities is not suitable for the analysis of the final states. In the next sections we dis- 
cuss the general factorization framework which makes use of the unintegrated parton correlation 
functions. All the details, including the derivation of the factorization and the definitions of the 
unintegrated parton correlation functions can be found in [3]. 



2 Integrated, unintegrated, fully unintegrated parton correlation functions 

We have learned from the discussion in the introduction that the description of the exclusive 
processes requires a new formalism with new parton densities and fragmentation functions. We 
distinguish three types of objects: 

• Conventional integrated parton distributions and fragmentation functions. These depend 
on the fraction of the longitudinal momentum and the scale 

:r(a;,/i) 

• Unintegrated parton distribution functions (and FFs) which depend on the transverse mo- 
mentum 

F{x, kt, fi) 

• Fully unintegrated parton correlation functions which depend on all components of the 
momenta: longitudinal, transverse, and virtuality 

The unintegrated (over kx) parton distribution functions have been introduced in the con- 
text of small X physics, (see for example [4]). There, as a result of the high-energy approximation, 
the 4-point gluon Green's function emerges as a fundamental object. It depends on the trans- 
verse momenta of 4-off shell gluons. Its evolution with rapidity is then governed by the BFKL 
equation. There have been numerous efforts [5-7] to construct the unintegrated parton densities 
outside the small x limit, many of them using the angular ordered CCFM equation. However, 
theses approaches do not provide with the operator definitions of the unintegrated parton dis- 
tribution functions. Therefore, it is desirable to have a unique and consistent approach which 
demonstrates factorization using the properly defined unintegrated parton correlation functions. 
Then, hopefully, different approaches developed so far will be reproduced as a particular limit or 
approximation to a more general case. 




Fig. 2: Most general graph with one jet production in deep inelastic scattering. 



3 General strategy 

The general formalism for use in the case of the deep inelastic scattering introduced in detail 
in [3]. The main points are 

• Retain the exact kinematics in initial and final states. 

• Keep the explicit factors (bubbles) representing final states. 

• Retain on-shell matrix elements. 

• Define projections from exact to approximate momenta. 

• Construct definitions of the gauge invariant parton correlation factors and the soft factor. 

• Use eikonal lines and Ward identities to prove factorization. 

• Use non-light-like eikonal Wilson lines to regularize light-cone divergences in the uninte- 
grated parton distributions ( this introduces a cutoff in rapidity, similarly to what is done at 
small X ). 

One starts with the most general graph in the case of the deep inelastic scattering which is 
illustrated in Fig. [3] This graphs differs from the ordinary parton model due to several distinctive 
features. First of all, the final state quark has a jet subgraph. This is a minimal required assump- 
tion since we know that the quark has to hadronize. A second important point is that we need 
to allow for the non-perturbative (or soft) interactions between the outgoing jet and the target 
remnant. These are obviously necessary to neutralize the color. Therefore we need to allow for 
arbitrarily many gluon connections between the soft bubble and the target and jet subgraphs. So 
graph in Fig. [3] constitutes a necessary extension to a simple parton model. 

In order to prove the factorization one needs first to identify all the leading regions: 
coUinear to jet, collinear to the target and soft. On also needs to perform subtractions to en- 
sure that the smaller regions are appropriately suppressed. Then the Ward identities are used to 
disentangle the gluons between the different subgraphs (for example collinear to soft subgraphs). 
In particular one has to prove that the Ward identities work properly in the presence of the sub- 
tractions. After the application of Ward identities the gluons factorize into the Wilson lines in the 
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Fig. 3: Most general graph with one jet production in deep inelastic scattering in the factorized form. Double lines 
denote the eikonal lines. 

appropriate objects. The final factorization formula proven in [3] reads 



where 5, JT, J ^ H denote the soft factor, parton distribution, jet factor and hard scattering cross 
section respectively. P^^^ denotes the projection onto the apropriate structure function. All the 
non-perturbative factors depend on the all components of the 4-momenta ki. Note the appearance 
of the rapidity variables yi. These are necessary to suppress the contributions from the regions 
where the rapidities are very large (rapidity divergencies). It has been demonstrated that this 
can be regularized by changing the direction of the Wilson lines from the light-like directions 
to slightly non-light-like [3]. The parton distribution, fragmentation functions and soft factors 
acquire then additional parameters. Appropriate evolution equations give the prescription on the 
variation of these quantities with the rapidity. The final factorized graph is shown in Fig. [3] where 
by double lines we denote the Wilson lines. 

4 Gauge invariant definitions of PCFs 

The analysis of the single gluon attachment allows to construct the gauge invariant definitions of 
the parton correlation functions (PCFs) and the soft factor. As already discussed in a previous 
section critical to these are the directions of the Wilson lines. Soft gluons couple to the target 
jet, with its large plus component of momentum, and to the outgoing jet with its large minus 
component of momentum. This means that in coordinate space the soft factor is the vacuum 
expectation value of Wilson lines that are nearly light-like in the plus and minus directions [3] 
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Note that we have replaced in this expression the Ught-Uke Wilson lines with the non-light-like 
ones in the directions denoted by the vectors ut = (1, — e"^^^^ Qt), and nj = (— e"'^^-'', 1, Ot) 
where yr ^ 1 and \yj\ ^ 1, yj < 0. As mentioned above, the reason for their presence is 
that they provide cutoffs on rapidity divergences. The factors / are the transverse links at infinity 
which are needed to ensure strict gauge invariance [8]. With these included, the soft factor (|7]) is 
the expectation value of the closed Wilson loop. 

The target PCF is defined as the gauge invariant expectation value of two quark fields. Its 
definition reads [3] 
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where w is a reference point in space-time. The two quark fields are connected via three Wilson 
lines, two of them going into the direction Us = (— e^', e^^% Ot) which corresponds to the 
rapidity close to zero in the center-of-mass system. This expectation value is divided by the soft 
factor, in order to ensure the proper cancellation of the double counting contributions. Note that 
the soft factor has now nj replaced by the vector. This guarantees that the definition ([8]) 
gives a good approximation for the gluons close to the target region. Similar definition can be 
constructed for the jet factor with vector tit replaced by n j 
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jiw, argument like yj, y,-},ys, = 1 1 — j • 
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